Linking Entanglement and Quantum Phase Transitions via Density Functional Theory 
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Density functional theory (DFT) is shown to provide a novel conceptual and computational framework for 
entanglement in interacting many-body quantum systems. DFT can, in particular, shed light on the intriguing 
relationship between quantum phase transitions and entanglement. We use DFT concepts to express entangle- 
ment measures in terms of the first or second derivative of the ground state energy. We illustrate the versatility 
of the DFT approach via a variety of analytically solvable models. As a further application we discuss entangle- 
ment and quantum phase transitions in the case of mean field approximations for realistic models of many-body 
systems. 



I. INTRODUCTION 

Density functional theory (DFT) |Q1 0] is to date the most 
successful method for first principies calculations of the elec- 
trònic properties of sòlids. The key to its success is a trans- 
formation of the dependence of the properties of a system 
of interacting partides on their single particle potential, to a 
dependence on the ground state density, thereby facilitating 
useful approximations of the many-body interaction for first 
principies computations. One rather relevant phenomenon in 
many-body physics is the ocurrence of quantum phase transi- 
tions (QPTs), which consist in critical changes in the proper- 
ties of the ground state, driven purely by quantum fluctuations 
and effectively occurring at temperature T = [3]. QPTs 
are associated with level crossings, which usually lead to the 
presence of non-analyticities in the energy spectrum. Specif- 
ically, a first-order QPT (1QPT) is characterized by a finite 
discontinuity in the first derivative of the ground state energy. 
Similarly, a second-order QPT (2QPT) is characterized either 
by a finite discontinuity or divergence in the second deriva- 
tive of the ground state energy, assuming the first derivative is 
continuous. 

Many-body physics and, in particular, critical phenomena 
near QPTs, have recently been the subject of intense interest 
from the perspective of the theory of quantum information. 
A key observation is that, since entanglement describes cor- 
relations in a quantum system, its quantification may provide 
a convenient and precise description of a QPT. Indeed, en- 
tanglement has been found to exhibit scaling behavior near a 
critical point fHHll^]. Moreover, under well-delineated con- 
ditions and for distinguishable systems up to two-body inter- 
actions, a formal relationship between QPT and bipartite en- 
tanglement was recently established [71]. Here, we show that 
entanglement may be well specified and conveniently com- 
puted within DFT. In DFT, any entanglement measure is a 
function(al) of the expectation vàlues of the observables. This 
procedure introduces a direct connection between entangle- 
ment and the derivatives of the ground state energy of the 



quantum system with respect to the field coefficients, leading 
to a deep relationship between entanglement and QPT. 



II. GENERALIZED HOHENBERG-KOHN (HK) 
THEOREM AND ENTANGLEMENT 

Consider a quantum system described by a Hamiltonian 
composed of two parts 

H = H +H eKt = H + J2*A, (1) 

i 

where A/ is the "field coefficient" (control parameter) asso- 
ciated with a set of Hermitian operators {^4/}, e.g., an ob- 
servable relevant to driving a quantum phase transition. The 
index l can be discrete or continuous. The expectation vàl- 
ues of Ai for a ground state \ip) are denoted by the set 

{ai} = {(->P\ Ai 

DFT is originally based on the Hohenberg-Kohn (HK) the- 
orem OJJ] . In the case of a many-electron system, the HK the- 
orem establishes that the ground state electrònic density n(r), 
instead of the potential v(r), can be used as the fundamental 
variable to describe the physical properties of the system. In 
the case of a Hamiltonian given by Eq. ([TJ, the HK theorem 
can be generalized to the statement that there is a duality (in 
the sense of a Legendre transform) between the set of expec- 
tation vàlues {et/} (corresponding to n(r)) and the set of field 
parameters {A/} (corresponding to v(r)) |8fl. The commuta- 
tivity of the densities at distinct points, [n(r), n(r')] = for 
r ^ r', is a property of the original HK theorem. In a lattice 
system, we require that the physical observables {Ai} on dif- 
ferent sites are mutually commuting operators. This allows 
not only different observables on the same site, e.g. Sf, Sf, 
to be non-commutative, but also, for later use, endows a func- 
tion of observables on different sites with a single site locality, 
such as the set of two site operators, {AiAi +c ], l ranging over 
all sites and c being a constant. It follows from the Legendre 
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transform that the ground state expectation value of any ob- 
servable can be interchangeably viewed as a unique function 
of either {Xi} or {ai}. (See Appendix lAl for a simple proof 
of the HK theorem in a lattice). Such a general duality has al- 
lowed for the application of DFT in, e.g., interacting quantum 
spin systems [9]. Moreover, as we show below, it can provide 
a natural connection between entanglement and QPT. Indeed, 
using the Hellmann-Feynman theorem lldlll 111 . 

d W dM 

^ = m^^) = ^\A l \i>) = a l . (2) 

This means that the set of observables {dH/dXi} has a direct 
linear relation with {a;}. An example is the metallization of a 
semiconductorunderpressure to a value at which the band gap 
given by the discontinuity of the density functional derivative 
of the ground state energy goes to zero lfÏ2Í [Ï3ll . 

The HK theorem can be used to redefine entanglement mea- 
sures in terms ofnew physical quantities: expectation vàlues 
of observables, {ai}, instead of external control parameters, 
{A;}. Consider an arbitrary entanglement measure M for the 
ground state of Hamiltonian (Q}. We will focus here on bipar- 
tite entanglement, but our discussion applies equally well to 
multipartite measures. We then prové a central lemma, which 
very generally connects M and energy derivatives. 
Lemma. Any entanglement measure M can be expressed as 
a unique functional of the set of first derivatives of the ground 
state energy: 

dE 

M = M({ ai }) = M({— }), (3) 

assuming that the ground state is non-degenerate. 
Proof. Intuitively, the proof follows from the fact that, accord- 
ing to the generalized HK theorem, any ground state wave 
function l't) is a unique functional of {ai} and since 
provides a complete description of the state of the system, 
everything else is a unique functional of {ai} as well, in- 
cluding M, More formally, let us consider the case of pair- 
wise entanglement of qubits. The case of higher dimensional 
systems or multipartite entanglement is a direct generaliza- 
tion. Then: (A) My (entanglement measure between qubits 
i and j) is always a function / of the matrix elements of 
the 2-qubit reduced density matrix p.y : AEj = f(pij). (B) 
The matrix elements are combinations of correlation func- 
tions (afer 1 -) = Tr(af cr^p.y), where a,b = 0, ..,3, with 
tro — I (identity). This follows from an expansion of pij 
in the Pauli basis {afa^}. (C) From steps (A) and (B) it fol- 
lows that M = M((of oj)). However, by using the HK the- 
orem for non-degenerate ground states, any expectation value 
can be taken as a function of {a;}, since the wave function 
itself is a function of {ai} (see, e.g., Ref. [8]). Therefore, 
M = M{{aA), as required. ■ 

In Ref. [7], relations similar to Eq. (O, which connects en- 
tanglement and lQPTs, were established at the critical point 
for several examples of multi-particle systems, up to two-body 
interactions. In DFT, Eq. (0 holds for arbitrary systems, and 
not only close to the critical point. While the HK theorem 
is also applicable to degenerate ground states iTTsL [HÜl . not 



all linear combinations of densities corresponding to degen- 
erate ground states are permissible when implementing the 
variational principle 11711 . Note also that systems described 
by either Fermi or Pauli operators can be considered using 
DFT. Indeed, the treatment of both cases can be unified by the 
Jordan-Wigner transformation [14], with H, Ho, and H ext 
expressed in terms of linear combinations of generators of 
SU(2 N ), where N denotes the number of sites in the case 
of spins in a lattice, or the number of single modes for Jordan- 
Wigner fermions. 

Moreover, the HK theorem implies that one can split up the 
Hamiltonian (Q} in different ways. For example, a new Hq 
might include part of the sum J^^iAi- In our discussion, it is 
often convenient to focus on one of the external operators by 
moving the others into Ho- 

For 2QPTs, we should examine the derivatives of M. For 
simplicity of exposition, we regard one of the parameters A; as 
an independent variable, which we denote by A, and consider 
all the others as part of Ho. Therefore, M can be seen as an 
exclusive function of A, yielding via Eq. (O 

dM dM da dM d 2 E 

= = . (4) 

d\ da dX da d\ 2 

Notice that this equation holds only for non-degenerate 
ground states, since for the case of degeneracy, although the 
density a still uniquely specifies the potential À, the potential 
A does not uniquely specify the density a anymore. There- 
fore, in the degenerate case, a cannot be taken as a function 
of A, which implies that the chain rule used to take the deriva- 
tive in Eq. is not vàlid. However, as long we restrict our- 
selves to non-degenerate states (as is the usual case for large 
finite systems tending to criticality), or approach the (criti- 
cal) degeneracy point from below or above, this problem can 
be avoided. In the case where the degeneracy is symmetry 
driven, we could also circumvent this problem by observing 
that degenerate states can be split by a symmetry breaking 
term which is then allowed to tend to zero in the study of 
QPTs. 

Eq. (0]i shows that an entanglement measure is proportional 
to the second derivative of energy as long as dM /da ^ 0. 
By using appropriate bipartite entanglement measures, 2QPTs 
have usually been identified so far through either non-analytic 
or vanishing vàlues of dM/dX at the critical point. Both cases 
are contained in Eq. ©. 

It should be emphasized that Eqs. (O and hold for any 
system described by the Hamiltonian (fl~|i as long as DFT is 
vàlid, in the degeneracy sense discussed above. Around the 
critical points, the left and right limits of the two equations 
still hold even if the DFT is qüestionable at the critical point. 
Eqs. (0 and (0 can be seen as the bàsic equations for the 
relation between QPTs and entanglement. 

III. EXAMPLE I: ONE-BODY EXTERNAL COUPLINGS 

As a first example of the applicability of Eqs. (01 and ©, 
let us consider ií cx t = J2i ^ < ' which represents a sys- 
tem of qubits acted upon via independent single-qubit con- 
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trol terms. According to DFT, the energy is a functional of 
matrix elements of one-spin reduced density matrices, i.e., 
E = where ~$ i = l$\ lf t \tp) = ^ t E is the 

Bloch vector, with components pf = af We con- 
sider a bipartition of the system, splitting it up into two 
parts. Then, assuming that the system is in a pure state, we 
can use the linear entropy as a measure of block entangle- 
ment, which reads — [d/ {d — 1)](1 — Trp 2 ), where 

< < 1 and p denotes a d x cí-dimensional den- 
sity matrix. Explicit computation of the block entanglement 

(2) 

of one qubit (the iíh) with the rest of system yields L\ ' — 

1 — IT^J 2 = 1 — S7iE ■ \/iE, which is a function of the pa- 
rameters A . In the case of fermions, we replace the Pauli ma- 
trices by fermionic operators according to the Jordan-Wigner 
transformation. Then, L^ 2 * 1 = 1 — (dE/dX Z i) 2 (number con- 
servation law for fermions implies the vanishing of dE/d\ X i 
and dE/d\ y i). 

In the case of a 1QPT, characterized by a discontinuity in 

— y (2) 

V i E, we have a corresponding discontinuity in L\ ' unless 
SJ t E ■ Vi-E" is continuous. Therefore, in this case, when 
all A i are taken as independent external parameters, the en- 

(2) 

tanglement measure L\ ' is an analytic function of the first 
derivatives of the energy, yielding a natural relationship be- 

(2) 

tween lQPTs and L\ . A general discussion of 2QPTs is, on 
the other hand, not as straigthforward, since the structure of 

(2) 

the derivatives of L\ will depend on the details of the model. 
Thus, it turns out to be more useful to analyze a concrete ex- 
ample. Let us consider the transverse field Ising chain, where 
H = — ^^ 1 ((Tfcrf +1 +Acrf), with iV denoting the number of 
spins along the chain and with cyclic boundary conditions as- 
sumed. In this model, a discussion of entanglement as a func- 
tion of the coupling A was first presented in Refs. 14j,|5[]. Due 
to translational symmetry we have p z — a z \ip] 



where e = E/N . Therefore 



3L< 2 > 



di< 2 > d 2 



'ÚL· 



-. Divergence 



dX ~ dpz dX 2 

of d 2 e/d\ 2 at the quantum critical point A = 1 will thus re- 
sult in that of dL^ /dX unless dL^ /dp z = 0, which is not 
the case in this example. This is demonstrated in Fig.Q] where 
we plot as a function of p z . Both the maximum and the 
singularity of the derivative occur at the critical point. We 
can also apply the DFT approach to pairwise entanglement 
measures. For instance, let us consider entanglement between 
nearest-neighbor pairs in the transverse field Ising model as 
measured by the negativity Af (TU]. From Eq. we have 
§T = §jfm?- Notice that the divergence in d 2 e/dX 2 at the 
critical point naturally leads to a divergence in dN / dX, since 
dJV/dp z is a non-vanishing function at the QPT, as shown 
in Fig. H] In fact, the maximum of dAÍ/dp z approaches the 
critical point as the number of sites increases. 



IV. EXAMPLE II: TWO-BODY EXTERNAL COUPLINGS 

In the case of two-body external couplings, we take ií ox t — 
J2ij A^-cr^cr", where a S {x, y, z}. This Hamiltonian repre- 
sents a system of qubits controlled externally via two-body 




FIG. 1: Block entanglement i/ 2 ' as function of p z for the transverse 
field Ising model. The maximum occurs at the quantum critical point, 
where p z ~ 0.6366. 
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FIG. 2: Derivative of the negativity with respect to p 2 for the trans- 
verse field Ising model. The maximum approaches the quantum crit- 
ical point, where p z ~ 0.6366. 



interactions. Entanglement between qubits i,j and the rest 
of the system can then be computed by taking the linear en- 
tropy for d = 4. This yields = 1 - |E(^-) 2 - 

We now analyze the behavior of in some important mod- 
els exhibiting QPTs. For example, for the XXZ spin chain, 
we have H = (-1/2) £f (afaf +1 + + Aa?af +1 ), 

where cyclic boundary conditions are assumed. The external 
Hamiltonian is taken as ií C xt = — (A/2) ^ ofcf+i- Direct 
evaluation of L^ 4 ' then yields (See AppendixIBl 



£(4) 



A 2 

2 



de_ 
dA 



s 
~2 



eA 



de_ 
dA 



, (5) 



where e = E/N. Notice that L (4) is a function of the DFT 
variable a = (crfof +1 ) = —2{de/dA) since, due to the 
HK theorem, the energy density e can be taken as a func- 
tion of a. Thus, discontinuities in (de/dA) will be directly 
reflected in This model exhibits two distinct QPTs, 
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which occur at A = 1 and A = — 1. In order to eval- 
uate we consider the ground state wave-function with 
vanishing magnetization, which favors the presence of entan- 
glement in the system. A 1QPT occurs at A = 1, which 
separates a ferromagnetic phase from a gapless quasi-long- 
range-ordered phase. At this ferromagnetic critical point, the 
energy density as N — > oo is continuous, and is given by 
e(A = 1) = -1/2 |[Ï9l |20h. However, its first deriva- 
tive is discontinuous, with (9e/9A)a^i+ — » — 1/2 and 
(d£/dA)/^_,i- — > 0. From Eq. (|3), we can see that this 
discontinuity is immediately manifested in the entanglement 
measure, since jumps from 2/3 to 5/6 at A = 1. A 
continuous QPT in the XXZ chain occurs at A = — 1, sep- 
arating the gapless quasi-long-range-ordered phase from the 
antiferromagnetic phase. For this case, it is useful to com- 
pute the first derivative of with respect to A, which 

yields (dL^/dA) = (dL^ /da)(d 2 e/dA 2 ), with ^ = 



eA- 2 



Se 



The QPT in this case is not di- 
rectly signalled by (d 2 e/dA 2 ), which is analytic at A = 

-i du. 

However, entanglement detects this transition 
as an extremum at the critical point ll2Ïll22l l23ll . This behav- 
ior is also reflected in terms of the DFT variable a. We have 
e(A = -1) = 2 (In 2 - 1/4) |Ï9ll20h, and fïnd for the first 
derivative of the energy (9e/9A)a— ►- i ~ 0.2954. Therefore, 
we obtain (dL^/da) = (dL^/dA) = 0. 

We now analyze the behavior of in a Fermi system. 
An interesting example is then the one-dimensional Hub- 
bard model, ií cx t = UJ2 a n at n ai, where n a [ (n a {) is 
the spin down (up) electrònic number at site a. The Hub- 
bard model describes a metal-insulating transition, which has 
been considered from the point of view of entanglement in 
Refs. Í24l l25ll . We can rearrange the indices for the modes 
a t and a j into nearest neighbor indices í and i + 1, re- 
spectively, in a linear lattice. Therefore, the Hamiltonian can 
be written as H cxt = U J2i n ii-\ n 2i where only the pairs of 
sites (1,2), (3,4), etc, interact with each other. We can then 
compute between an interacti ng p air (i, i + 1) and the 
rest of the system (see also Refs. |24Í I25I0 . At half-filling, 
= §(1 + 4a - 8a 2 ) (for any i) (See Appendix 0. 

Then = |(4 - 16a). By using Eq. ©, we obtain 

= ^Ttrmz- At the critical P oint U = 0, which 
separates an insulating phase from a metallic phase, the first 
derivative of L {4 > with respect to U is dL^/dU = Q. 
In terms of the new variable a, we can show that the QPT in 
the Hubbard model is also identified via an extremum of . 
Indeed, for U = 0, we have a = 1/4 Í26ll which then implies 
d a LW = 0. 



V. THE LIPKIN MODEL : A HARTREE-FOCK APPROACH 
TO ENTANGLEMENT 

Most realistic physical many-body problems cannot be 
solved analytically. Linear approximations, such as Hartree- 
Fock-Bogoliubov theory and the random phase approxima- 
tion, are often practical and effective ways to treat these 
systems, since these procedures change an intractable 2^- 



dimensional problem to a tractable A f2 -dimensional one. In 
this case, it is appealing to introduce new and simple quanti- 
ties, e.g., and as measures characterizing the quan- 
tum information content of these known approximate wave 
functions. We expect these quantities to become as impor- 
tant as, e.g., binding energies, when quantum information be- 
comes readily accessible to experiments. 

As an example, we consider the Lipkin model - impor- 
tant, e.g., in nuclear physics - whose Hamiltonian reads 
H = XS Z — jj(S 2 — S 2 ), where S z = J2m=i è( c +™ c + m — 
cL m c_ m ) and S x + iS y = X)m=i c +m c -m (for a 
discussion of entanglement in the Lipkin model, see also 
Ref. lEÜO . This Hamiltonian describes a two-level Fermi sys- 
tem {|+), |— )}, each level having degeneracy N. The oper- 
ators c , m and cL m create a particle in the upper and lower 
levels, respectively. Alternatively, the Hamiltonian may be 
viewed as a one-dimensional ring of two-level àtoms with 
infinite range interaction between pairs. The factor 1/N in 
the interaction term keeps the scaling of both terms in H 
linear in N. The phase transition studied is in the limit 
of N — * 00. The Lipkin model is exactly solvable (see, 
e.g., Ref. 1I29I0 . The Hartree-Fock (HF) ground state, which 
is exact for this model as N tends to infinity, is given by 
\HF) = Y[<m=i a oml~)' wnere a lm is defined by the follow- 



ing change of variables: 



and 



- +m - smaa„ m + cosaai m 
c_ m — cos aal m — sinaa\ m . The variational parameter a 
which yields the minimum energy is given by cos 2a = X 
when A < 1 and a = when A > 1. We define the 
DFT variable a = de/dX, with e — E/N denoting the en- 
ergy per particle. For the HF ground state, we then obtain 



for A < 1 and a 



for A > 1 It is easy to 



show that d 2 e/dX 2 is discontinuous at À = 1, which corre- 
sponds to a = —1/2 in terms of the DFT variable. Let us 
analyze whether this discontinuity is reflected in the deriva- 
tives of the entanglement measures, as given by Eq. For 
4-dimensional block entanglement, it is convenient to con- 
sider the entanglement between a block composed of two 
general modes (+m,— n) and the rest of the system, which 
yields Lf m _ n = (2/3)(l - 4a 2 )(l - 6 m , n ), where 6 m , n 
is the Kronecker symbol. Therefore, the block (+m, —n) 
is entangled with the rest of system only if m / n. Tak- 
ing the derivative, we obtain (9L^ _„/9a) a= _ 1 / 2 = 8/3 
(m 7^ n). Therefore, from Eq. (0}, the non-analyticity of 
d 2 e/dX 2 at the critical point will be associated to a non- 
analyticity in _„/9A (m 7^ n). A similar result fol- 
lows in the case of 2-dimensional block entanglement, where 
we have = 1 — 4a 2 for a general mode +m (or — m) 
with the rest of the system. Pairwise entanglement between 
general modes +m and — n as measured by the negativity is 
found to be A/+ m ,-n = Vl — 4a 2 S m , n - Notice that this is in 
contrast with block entanglement, where modes +m and — n 
only are entangled for m ^ n. This difference is due to the 
structure of the HF ground state, which implies that the modes 
+m and — n interact only for m = n. Therefore, bipartite 
entanglement in the system appears only when +m and — m 
are in different parts. Evaluating now the derivative of the 



5 



negativity we obtain (dN +m ,- n / da) a ^_ l/2 -> 
dAÍ+m^-n/dX is non-analytic at the critical point. 



VI. CONCLUSION 



oo. Thus, with a[ = (ip' \Ai\ip') . From Eqs. iA3\ and (IA41 > we have 

< E ( A í - A; ) " °í) (A5) 



We have shown in general and illustrated in a number of 
models that DFT provides a natural link between entangle- 
ment and QPTs. Since experimental data are taken at finite 
temperature, it is important to be able to delineate the tem- 
perature fluctuation around a classical critical point versus the 
quantum fluctuations around a QPT. The exploration of finite- 
temperature DFT 13011 for the connection between phase tran- 
sitions and quantum information appears to be a promising 
direction for future study. 
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APPENDIX A 

We provide here a proof of the HK theorem in a lattice, 
based on the variational method (for a proof based on the 
constrained-search technique B15I1 . see Ref. (8|]). Let us con- 
sider two sets of parameters {A/} and {AJ}, which define two 
Hamiltonians as follows: 



H = H 



A/A; , H = H 



(Al) 



The ground states of H and H' will be denoted by \tp) 
and \ip'), respectively, which are taken as non-degenerate, 
even though the proof can be extended for degenerate ground 
states BlSLIlol . We also assume here that, for different sets of 
parameters, {A;} ^ {AJ}, we have independent ground states 
\tp) ^ a\ip') (a = constant). This is indeed the usual behav- 
ior of quantum systems around criticality, where the ground 
state varies continuously as we vary the control parameters. 
By applying the variational principle for the Hamiltonian H' , 
we obtain 

(v' W) < (m'w = a>\ [ H+ T, ( A í - A <) w 

(A2) 

Therefore, Eq. ( IA2t yields 

E' <E + J2M-*i)ah (A3) 
i 

where E and Eq are the ground state energies of H' and H, 
respectively, and a\ — (ip\Ai \ip). Analogously, by applying 
the variational principle for H, we obtain 



EoKE'c + ^iXt-X'Aal, 



Hence, if the sets of parameters {A;} and {X[} are different 
from each other, then we cannot have identical sets {ai} and 
{a[}. Therefore, the density {a/} uniquely specifies the poten- 
tial {A;} and can then be used as the bàsic variable to describe 
the properties of the system. 



APPENDIX B 

We provide here the bàsic details of the evaluation of the 
linear entropy for the XXZ model. The density matrix for a 
pair of nearest-neighbor sites in the ground state of the XXZ 
chain can be written as 



P = 



A 
OSCO 
C B 
-D 



(Bl) 



where, from the XXZ Hamiltonian, we obtain 



1 / de 

A = D = -;\ 

4 V <9A 



C 



1 ( A Qe 

"2 [ £ - A dA 



(B2) 



Eqs. (IB2t allows for a direct calculation of the linear entropy 
£( 4 ) = (4/3)(l — Trp 2 ), yielding the result displayed in 
Eq. ©. 



APPENDIX C 

We provide here the bàsic details for the evaluation of 
the linear entropy in the Hubbard model. Translation in- 
variance and simultaneous conservation of particle number 
N = J2j ( n j1 + n jl) anc l -z-component of total spin S z — 
J2j (nj-\ — riji) imply that the density operator for any sin- 
gle site can be represented by a 4 x 4 diagonal matrix, whose 
eigenvalues are given by 



vT = (n al ) 



de_ 
dÜ 

- w, z = 1 — u + — u~ — w. 



a, u — {n a f) — w, 



(Cl) 



(A4) 



At half-filling, we have (n a f) — (n a j_) = 1/2. Therefore, in 
this regime, all the eigenvalues can be expressed in terms of 
the density a = (n a ^n a {). Then, the evaluation of the linear 
entropy L^ 4 \a) follows straightforwardly. 
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